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Abstract 

The most general black M5-brane solution of eleven-dimensional supergravity (with 
o ■ a flat R^ spacetime in the brane and a regular horizon) is characterized by charge, 

O . mass and two angular momenta. We use this metric to construct general dual models 

of large-A^ QCD (at strong coupling) that depend on two free parameters. The mass 
spectrum of scalar particles is determined analytically (in the WKB approximation) 
^ ■ and numerically in the whole two-dimensional parameter space. We compare the mass 

spectrum with analogous results from lattice calculations, and find that the super- 
gravity predictions are close to the lattice results everywhere on the two dimensional 
parameter space except along a special line. We also examine the mass spectrum of 
the supergravity Kaluza-Klein (KK) modes and find that the KK modes along the 
compact D-brane coordinate decouple from the spectrum for large angular momenta. 
There are however KK modes charged under a U{1) x U{1) global symmetry which 
do not decouple anywhere on the parameter space. General formulas for the string 
tension and action are also given. 
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1 Introduction 

The conjectured dualities between gauge and string theories [1] have been recently exploited 
in [2]- [9] to construct and investigate models of pure QCD in 3+1 dimensions, whose main 
component is the black M5-brane solution of eleven-dimensional supergravity, which near the 
branes corresponds to an anti-de Sitter (AdS) space. The no-hair theorem implies that the 
most general model of this kind that can be constructed (i.e. based on a regular geometry 
with M5-brane charge) is obtained from a rotating black M5-brane parameterized by its 
charge, mass and two angular momenta. The scope of this paper is to calculate the mass 
spectrum of scalar modes of this general model in the supergravity approximation, and study 
its behavior in the parameter space. The parameter space is four dimensional, but the mass 
parameter can be set to 1 by a choice of mass units; the charge is related to the 't Hooft 
coupling A = g'^N (where g is the Yang-Mills coupling and A^ is the number of branes). 
It is assumed that A is very large so that the radius of curvature is much larger than the 
string scale; this is necessary for supergravity to be a good approximation to string theory 
(M theory). In this regime glueball masses are independent of A, so what remains is a two- 
dimensional space spanned by the angular momentum parameters. When one of the angular 
momenta vanishes, the model reduces to the one angular momentum model examined in 
Refs. [6, 8]. In our investigation we will use both analytic methods (within the WKB 
approximation, as developed in Refs. [9, 10]) as well as numerical ones based on Ref. [4]. 

The static M5-brane has an 5*0(5) symmetry associated with the internal S*^. Turning 
on the angular momentum parameters, this symmetry group breaks down to the Cartan 
subgroup as 5*0(5) — > 5*0(2) x 5*0(2) ~ U{1) x U{\). The spectrum of the supergravity field 
fiuctuations can be organized in representations of 5*0(5) or 5*0(2) x 5*0(2). The proposal of 
Refs. [2]-[4] is to identify the 5*0(5)-singlet modes propagating on the Minkowski boundary 
of the spacetime with large- A^ QCD glueballs. The dilaton modes correspond to J^'-' = 0'^"'" 
glueballs (J, P, and C being the spin, parity and charge conjugation quantum numbers). In 
non-supersymmetric, pure SU{N) Yang-Mills theory, there is no counterpart of the 5*0(5) 
global symmetry, so one would expect that at weak Yang-Mills coupling those Kaluza-Klein 
(KK) particles which transform non-trivially under this group are very massive and decouple. 
This problem was studied for QCD3 in [11] where it was shown that the first correction 
(beyond the A = 00 limit) to the masses of these states does not lead to their decoupling in 
the case of vanishing angular momenta. A general study for QCD3 supergravity models with 
three angular momenta was recently given in [10]. In this paper, using both analytic (within 
the WKB approximation) and numerical methods, we calculate the spectrum of glueballs 
and of KK states. We find that the KK modes on 5*'* do not decouple in the large A regime in 
any region of the two dimensional parameter space (within the supergravity approximation). 
In contrast, the KK modes on the circle associated with the compact Euclidean time (on the 
M5-brane worldvolume) decouple in the limit of large angular momentum. 

Some interesting effects concerning thermodynamical aspects of rotating D-branes have 



been recently pointed out in Refs. [12]-[14]. Here we will be considering the slightly different 
construction of Refs. [2, 6] for zero-temperature QCD, where the Euclidean time parame- 
terizes an internal circle, and the Minkowski time is one of the brane volume coordinates. 



2 The Supergravity Model 

2.1 The metric 

The maximal number of angular momentum parameters for the rotating M5-brane (dictated 
by the rank of the 5*0(5) isometry group of rotations of the static M5-brane) is equal to two. 
This metric was constructed in [15], though the expression given there contains a few minor 
mistakes which we correct below. The metric of the rotating M5-brane is given by^ 

,2 



dsl, = rH-hdt' + dxl + ... + dxl) + ^dr^ + fir%l + '-^^^ + 



2x /^ 7 2 r^ 2r/-, ll COS^ 9 

h r^ 



-)d9^ + (1 + - — -— -)cos^6'rf^^-2^cos6'sin6'cos^sin^ci6'# 



/jsin^^sin^V;^^^, ^ ^^ ^ Ijcop^^ ^^^, ^^^, _ ^l? 

cosh5 ,,, .2/17 1 2/1-2,7 \ 4m/i/2 cos^6'sin^6'sin^-?/' 
~4m— r— — -atf /i sm udipi + L2C0S dsvn ipdip2) -\ rin- dipidip2 

o ^ , /? 2m/?sin^^, , , ,^ , , , /, 2772/, cos^ ^ sin^ -^ , , 
+ sm2^(l + -1 + ^l^^ )d^l + cos2^smV(l + -| + '—^. ^)rf^2j 

(2.1: 



where 



/2 72 72/2 

A = 1 + -^ cos^ B + %(sin2 e + cos^ 9 cos^ ^) + ^ cos^ 9 cos^ ^ , (2.2) 





cos^ 6' cos^ 


^)- 


^ ^ ^ 2m sinh^ a 
^ - 1 + Ar3 


5 


/i=l 


2m 
Ar3 ' 




_ 1 + ^ + 


,2 ,2,2 

h _|_ 'i'2 


2m 
^3 



(2.3) 
(2.4) 

The horizon is located at r = th-, where th is the largest real root of 

(r2 + /^)(r^ + /^)-2mr = . (2.6) 

One can obtain the following formulas for the ADM mass, entropy and angular momentum: 

^This differs from Eq. (12) of [15] in the expression for A (called /^^ there), the power of r in the 
components gtipi,gtip2, 3viV2: a-nd a factor sin^ ip in gt^^. 
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where G^ is Newton's constant in 11 dimensions, and /p is the 11 dimensional Planck length. 
The parameter a is related to the (magnetic) charge N and m by 



sinh^ a = -{ ^{-KNll/mf + 1-1 



The Hawking temperature and angular velocities are given by 
_ 3r|, + {ll + llYl - llll 

^H — 5 '-, 1 "1,2 



STfmrfj cosh a 



cosh a {rjj + If 2) 

These quantities satisfy the first law of black hole thermodynamics: 
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Let us now go to Euclidean space r = —it, Zi 2 — * '^/i,2j and take the field theory limit as 
in [1, 6]: 

r = UHl, 2m = U^ll, /i,2 = ajjl , (2.14) 
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af cos^ ^ a|(sin^ ^ + cos^ ^ cos^ '?/') a\a 2 coip' 6 cos^ ip 
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Note that the component gLp-^Lp2 vanishes in the field theory hmit, and so do the last terms 
^v^iv^i <ino. yip2^2' 

The coordinate r describes a circle of radius Rq, where Rq is related to the Hawking 
temperature Th by Rq = {2ttTh)~^, with 
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where we have introduced the coordinate u by U = 2{tx NY^'^u^ and rescaled 01^2 -^ 
2(7riV)^/^ai 2. The constants u'jj,ujfj,u1,U2 represent the four different solutions for -u^ of 
the equation 

{u" - at) {u" - 4) - uy = . (2.21) 

There are two positive {ujj, u]^, ujj > uj^j-), and two negative (or complex) solutions (m^, u^), 
with Uff and ujj^ representing the outer and inner horizons respectively. When ai = 02 = a, 
the equation simplifies to 

M^ - a^ = ±ulu , (2.22) 

where the signs ± corresponds to the inner and outer horizons. From Eq. (2.22) one sees 
that when a ^ Uq the two positive solutions get closer to each other, thus the inner horizon 
approaches the outer horizon. 

The gauge coupling gl in the 3 + 1 dimensional Yang-Mills theory is given by the ratio 
between the periods of the eleven-dimensional coordinates X5 and r, i.e. 

r = i?o^2, X5 = |^i?o^i = 4^o^i . ^1,2 = ^1,2 + 27r , (2.23) 

where A = ^^ is the 't Hooft coupling. Dimensional reduction in 61 gives the type IIA 
metric representing the field theory limit of the rotating D4-brane metric with two angular 
momentum parameters: 

^^^ 3Uq L V U i 2 d/ g^2 ^6^^ 2. 
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where the dilaton field is given by 

In these coordinates, the metric is independent of A^, and the string couphng is of order 1/A^, 
as expected. The 't Hooft couphng A appears as an overall factor of the metric. For uq 7^ 0, 
curvature invariants have a finite value at the horizon, and they are suppressed by inverse 
powers of A. 

The metric (2.24) has a f/(l)'^ isometry associated with translations in 62,ipi,^p2- This 
should appear as a global symmetry in the corresponding dual Yang-Mills theory. Since the 
pure SU{N) QCD has no such symmetries, one may expect that states which have charges 
with respect to U{1)^ have a large mass compared to the glueball masses. In Section 4 we 
calculate the different mass spectra and investigate this possibility. 

2.2 String tension and action 

The string tension is given by l/27r times the coefficient of J^dxf, evaluated at the horizon, 
at the angles where it takes its minimum value [2, 6]. This follows by minimizing the Nambu- 
Goto action of the string configuration. The absolute minimum occurs at ^ = ■(/' = or tt. 
We obtain 

a = -XAul . (2.26) 

o 

String excitations should have masses of order cr^/^. The spin < 2 glueballs that remain in 
the supergravity approximation - whose masses are determined from the Laplace equation - 
have masses which are independent of A. 

In the field theory limit, the free energy F (= Action x Th) takes the simple form 

F = E-THS-niJi-n2J2 = -^ N^ul, (2.27) 

where E = Madm - Mextremai , Mextremai = Madm(mo = 0). Usiug that the M5-brane coordi- 
nate Xs is compactified on a circle with radius RqX/N, one has the relation 

^5 = 7^ . (2.28) 



Expressing -Uq in terms of the string tension (2.26) we obtain the intriguing relation 

-a'^ , (2.29) 



Action 1 N"^ 2 
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that generahzes the result found in [8] for the case of one angular momentum. Thus, in terms 
of the string tension, the action is independent of ai^2- It would be very interesting to have 
a derivation of (2.29) from the Yang-Mills side as a non-perturbative contribution to the 
partition function (related to the expectation value of the gluon condensate (j-j — Tr -Ffj,(0))). 

3ym ^ 

2.3 The supersymmetric limit wq = 

Metrics of rotating branes with non-extremality parameter m = greatly simplify upon 
introducing Cartesian-type coordinates [6]. For the extremal {m = 0) M5-brane metric 
(2.1), one introduces [13] 



1/1 = Y r^ + /f sin 6 cos (fi , 1)2 = y ^^ + li sin 6 sin Lpi , 

yz = Y ^^ + ^2 cos 6 sin if) cos (/92 , ^4 = Y ^^ + ^2 cos 6 sin ■?/' sin ^92 , 

yr, = rcosOcosip . (2.30) 

Using these coordinates we obtain 

ds'uA = r"'[ -dxl + jz d^l] + f" E dy] , (2.31) 

i=i i=i 

where / is obtained from Eq. (2.3) by taking the limit a -^ oo, m ^ at fixed A^ using 

(2.11): / = 1 + -^^, with r, ^, ip expressed in terms of yj by Eq. (2.30). In this limit 
the BPS bound is saturated, Madm = const A^. It can be shown that the function /(|/j) 
satisfies the equation djd^ f = 0, i.e. it is a harmonic function in the 5-space parameterized 
by yj. The metric (2.31) has unbroken supersymmetries, which can also be understood by 
interpreting the metric (2.31) as a multicenter distribution of BPS D4-branes, by constructing 
the harmonic function / as a linear superposition of harmonic functions corresponding to 
each D4-brane [13, 16]. 

The field theory limit of (2.31) can be obtained by replacing / — > / — 1, and properly 
rescaling coordinates. Alternatively, we can return to the metric (2.24) written in spherical 
coordinates, and set Uq = 0. The resulting metric has a curvature singularity in -u = ai (we 
are assuming ai > a2 > 0), which cannot be removed by any choice of periodicity in the r 
coordinates (the horizon region of the extremal uq = metric is not a Rindler space). Because 
of the singularity, the supergravity approximation breaks down in the uq = case; in order 
to understand the corresponding supersymmetric gauge theory, one needs to understand the 

6 



full string theory. At the supergravity level, it is meaningless to associate a temperature to 
this metric. 

One can have control over the string-theory corrections if we regularize the metric by 
taking mq 7^ and consider the limit of small uq (or equivalently, ai^2/uo large). For any value 
of ai^2/uo, one can choose A sufficiently large so that all curvature invariants are arbitrarily 
small. This is the technique used in the next section when discussing the large 01^2/^0 limit. 
Note that in this limit T^ ^ 00. In this theory, all fermions — which have masses 0{Th) — 
decouple. On the other hand, the spectrum of the Uq = theory must be supersymmetric 
with the usual degeneracy between fermions and bosons. The supersymmetric uq = theory 
cannot coincide with the theory obtained by taking the limit uq —^ (in the fashion described 
above), since the latter theory does not have fermions in the spectrum. Nevertheless, since 
the corresponding background metrics are essentially the same (at 01^2 ^ ""o) it is possible 
that a part of the structure of the theory with ai/uo ^ 1 may be dictated by the structure 
of the Uq = supersymmetric model. We shall return to this point in our conclusions. 

3 Glueballs and the Related KK Modes 

The 0"*"^ glueballs are related to spherically symmetric modes of the dilaton fluctuations, of 
the form 

^ = (l){u)e"'-^ , (3.1) 

where M^ = — /c^ [2]. The differential equation determining the mass eigenvalues is obtained 
by substituting this into the dilaton equation of motion 



V9 



e'^'^g^'^d^^l = , (3.2) 

using the background metric (2.24), and the formula 

v/^ = CM'^Acos2^sin^sin^ , C = -^ 1^^^] . (3.3) 



27rA 
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In addition to the 0"'""'' glueballs we consider particles with non- vanishing f/(l) charge asso- 
ciated with the circle parameterized by 62- The corresponding solutions of (3.2) will be of 
the form 

^ = 0(M)e*^-^e™''2 . (3.4) 

We will show both analytically (within the WKB approximation) and numerically that these 
states do decouple for a particular range of parameters. We will also consider the KK states 
associated with the / = 1 modes of the S"^. For the static (ai = 02 = 0) M5 metric, these 
transform in the 5 representation of 5*0(5). After introducing angular momentum, this 
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decomposes into (2, 1) © (1, 2) © (1, 1) of the Cartan subgroup 5*0(2) x 5*0(2). According 
to this decomposition, the corresponding solutions of (3.2) for the two doublets will be given 
by 



^ = (f)(u)e"'-''sme{ r^' , (3.5) 

^ sm (/9i ' 



^ = 0(M)e*'=- cos ^ sin V^ "r ^' , (3.6) 

V sm (^2 / 

whereas for the singlet it is of the form 

^ = 0(M)e*^cos^cosV' . (3.7) 

In ordinary (finite A, A^) Yang-Mills theory there is no 5*0(2) x 5*0(2) symmetry, so one 
would expect that at least the states which transform non-trivially under 5*0(2) x 5*0(2) 
become very massive and decouple in the weak-coupling limit. It is clear that the singlet state 
(3.7) should also decouple. If it did not decouple at small A, it would then be represented by 
some (gluon field) operator in the gauge theory. In the zero angular momentum case, this 
state combines with the other four components to form a multiplet (a 5) of 5*0(5). Thus 
the singlet state cannot correspond to a purely gluonic operator (since the gluon field is a 
singlet under 5*0(5)), and must decouple. 

Finally, we shall also consider ^ glueballs, which couple to the operator O4 = Tr FF. 
On the D4-brane worldvolume, the field that couples to this operator is the R-R 1-form A^, 
which satisfies the equation of motion 

d,[y^g^''g''^id,A^ - d^Ap)] = , /x, z/ = 1, . . . , 10 . (3.8) 

Finding angular-independent solutions is complicated, because of the non-diagonal compo- 
nents of the metric. The metric becomes diagonal in the two opposite limits ai,2 <^ uq and 
^1,2 3> Uq. In these cases one can consider solutions of the form 

Ae, =xe.(«)e*% A^ = if fi^e,. (3.9) 

In the following we will first present the mass spectra of these states obtained in the 
WKB approximation, and then the same spectra obtained by using numerical methods. We 
present tables for each state comparing the WKB with the numerical results and find that 
they are in a very good agreement. We also compare them to the lattice results for the 
glueball states which were computed for A^ = 3 and small A. 

3.1 Mass spectrum in the WKB approximation 

In the following we use the WKB approach of [10] (which generalizes the WKB approach 
of [9]) to calculate the different mass spectra (including KK modes) in the present case of 



QCD4 with two angular momenta. Consider differential equations of the form 

dp {f{p)dp(t>) + {m^Kp) + p(p)) = 0, (3.10) 

where M represents a mass parameter, and /(p), h{p) and p(p), p G [ph.oo), are three 
arbitrary functions which are independent of M and have the following behavior: 

I^Mp-PhY', h^hi{p-pHy\ P^Pi{p-pHy\ as p^pH, (3.11) 

/ ~ /2P''' , h^h2p''\ p^p2p''\ as p^oo, (3.12) 

where ^1^2,3, ■51,2,3, /i,2, ^1,2 and pi,2 are (real) numerical constants. For large masses M, the 
WKB method can be applied to obtain the approximate spectrum. One finds [10] 

TT^ / / ^2 f3- 

—rm{m+[-l^ h — 

where 



M2 = ^mfm+(-l + — + ^) I + C(m°) , m>l. (3.13) 



e = / ^p J7 , (3.14) 

is a constant which scales like a length, and 

ai = S2 - si + 2 , /3i = ri - r2 - 2 , (3.15) 

I Pi 

a2 = |si-l| or ^2 = W(si - 1)^ -4— (if S3 - Si + 2 = 0) , 



/52 = |ri-l| or /52 = ^(ri-l)2-4^ (ifri-r3-2 = 0) . (3.16) 

Consistency requires that ai and Pi are strictly positive numbers whereas S3 — si + 2 and 
^1—^3 — 2 can be either positive or zero. Typically the validity of the WKB approximation 
requires that the quantum number m be much larger than 1 (for precise conditions see [10]). 

3.1.1 Masses of the 0++ glueballs 

The masses of the O"*"^ glueballs are determined from the differential equation (3.2) with the 
ansatz (3.1). Introducing p = u^ one gets Eq. (3.10) with^ 



f{p) = (P' - &?)(P' - ^2) - PoP = {P- Ph){p - Pi)(p - P2)(p - P3) , 

Kp) = ^ , Pip) = . (3.17) 



'in the rest of subsection 3.1 we will use the notation af 2 = ^1.2 



This gives for the various constants 



Sl = 1 , 


S2 = 0, 


n =4 , 


r2 = l, 


ai = 1 , 


a2 = , 


/9i = l, 


/92 = 3. 



(3.18) 



Using (3.13) one obtains the following mass spectrum 
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M^ = — m(m + 2) + C(m°) , m > 1 , 
1 roo rfp ^ 



C=-/ '^ ^'^ (3.19) 

2 JpH y'(p - ph) (p - pi) (p - P2) (p - ps) 

This formula implies that mass ratios between resonances are, in the WKB approximation, 
independent of the angular momentum parameters &i,62- This is similar to QCD3, where 
the general rotating D3-brane solution with three angular momenta parameters was used 
[10]. As in [10], the WKB approximation breaks down in the region near bi = 62, Uq = 0. 

3.1.2 The KK modes on the circle 

For the KK modes with non- vanishing f/(l) charge corresponding to the periodic variable 62 
we look for solutions of (3.2) with the ansatz (3.4). In this case we obtain Eq. (3.10) with 
M^ replaced by M^ - Atc'^u'^T^ and 

f{p) = {p'-bl){p'-bl)-plp, 

Kp) = 4 , Pip) = 77-^ • (3-20) 

This gives for the various constants 

Sl = 1 , S2 = , S3 = -1 , ri = 4 , r2 = 1 , rs = -2 , 

ai = 1 , a2 = 7 -J 77 r , Pi = 1 , P2 = ^ ■ (3.21) 

[PH - pl){pH - p2){pH - P3) 

Using (2.19) and (2.20) we see that 02 = n. Then (3.13) (with M^ -^ M^ - Att^u^T^) gives 
the following mass spectrum 

M^ = An^n^Tl + — m{m + 2 + n) + 0{m^) , m > 1 . (3.22) 

We would like to examine the way that these states decouple in two limiting cases. First 
consider the case with bi ^ po and 62- Then we see from (2.19) that 

bl 



^ — 3/2 • 
Trpo 
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(3.23) 



On the other hand, in the same hmit, 

e - ^ • (3.24) 

Therefore the ratios of the masses of the glueballs to those of the f/(l) charged particles 
behave as 



^glucb. 



n \biJ 



.. , , , , for 61 > po and 62 • (3.25) 

Mcirc. n \hiJ 

Hence, the KK modes on the circle decouple with a power law. Now consider the case 

61 = 462 > po. Then 

156? 
Th ^ —4/2 (3-26) 

levrpo' 

and 

« =. i§ . (3.27) 

Therefore the ratios of the masses of the glueballs to those of the f/(l) charged particles 
behave as 

— ^ ^ 1.26-!^ — , for 61 = 462 > Po , (3.28) 

showing that in this case there is also decoupling with the same power law as in (3.25) (up 
to a slightly different numerical factor). 

3.1.3 The KK modes of S^ 

Let us now compute the mass spectrum for the KK states with non-trivial angular de- 
pendence on S^. The two equations corresponding to the doublets must be related by an 
interchange of Oi and 02, whereas the one corresponding to the singlet should be invariant 
under such an interchange. For the (2, 1) doublet we make the ansatz (3.5). Inserting into 
Eq. (3.2) one finds equation (3.10) with 

/(p) = (p'-&?)(p'-&2)-p^P, Kp) - P 



4 ' 
2pV f(p) V P\ 



,„)._vri_j|i_f^r,_|v. ,3.29) 



For the various constants we find 

si = 1 , S2 = , S3 = -1 , ri = 4 , r2 = 1 , r^ 
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fl = {pH - pl){pH - p2){pH - P3) , /2 = 1 , Pl = ij^ll- 



-'2 



h V pI) ' 

P2 = -4 , ai = 1 , a2 = 2^ [ 1 - -^ ) , /5i = 1 , /^s = 5 . (3.30) 

Ji V PhJ 

Hence, the mass is given by 

M2 = ^mL + 4+ ^(1-4)1 +C)K), rn>l, (3.31) 

where ^ is given by (3.19). Consider the mass formula in the region hi ^ po where the KK 
modes on the circle decouple, as in the case of one angular momentum [6, 8]. Using p^ = bi, 
the mass formula takes the form 

M^ ~ — m{m + 5) + C(m°) , m > 1 . (3.32) 

4 

This shows that for bi ^ po the mass of these KK states is of the same order as the glueball 
masses (3.19). 

For the (1, 2) doublet we make the ansatz (3.6) We obtain the same results as (3.29)- 
(3.31) with bi and 62 interchanged. For completeness we include the mass formula 

M2 = ^mL + 4+^(l-4)) +0{m'), m>l. (3.33) 

For ph — bi this becomes 

M^ ~ — m{m + 4) + C(m°) , m > 1 , (3.34) 

4 

where C, is given by (3.19). This shows that for 61 ^ po the mass of these KK states is 
of the same order as glueball masses and a little lighter than the modes corresponding to 
the KK doublet (3.5). As a general rule (which applies in particular to QCD3 [10]), states 
with (y9-dependence corresponding to the largest angular parameter are slightly heavier. In 
addition to the two doublets there is also a singlet (1, 1), represented by Eq. (3.7). We find 
that the function 0(p) obeys (3.10) with 

f{u) = {p'-bl){p'-bl)~plp, %) = ^, 

p(p) = 2{bl + bl- 2p2) . (3.35) 

For the various constants necessary to compute the corresponding masses we find 

Si = 1 , S2 = , S3 = , ri = 4 , r2 = 1 , ra = 2 , 

P2 = -4 , /2 = 1 , ai = 1 , a2 = , A = 1 , 132 = 5. (3.36) 
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Using (3.13) the mass formula for the singlet (3.7) reads 



vr^ 



M^ = — m{m + 4) + 0{m^) , m > 1 , (3.37) 

where ^ is again given in (3.19). Clearly, the masses of these modes are of the same order as 
the glueball masses (3.19), albeit slightly heavier. 

3.1.4 Masses of the 0+ glueballs 

Let us finally also consider ^ glueballs. As we have mentioned, finding angular-independent 
solutions is complicated, because of the non-diagonal components of the metric. The metric 
becomes diagonal in the two opposite limits 61^2 *^ Po ^-^d 61 2 ^ Po- Iii these cases one can 
consider solutions of the form (3.9). Substituting this into (3.8), we obtain a second order 
ordinary differential equation which, upon introducing p = u^ and writing Oi 2 = ^1,2 > can be 
written as Eq. (3.10) (with 0(p) -^ Xe2{p) ) with 



(3.38) 







fip) = ip'-bl)ip'-bl), 
,.. 1 p{p'-bl){p'- 
"^P^-Aip^-hDip^-hl) 


P{p) - 
hi) 

- Pop ' 


= 


This gives for the 


1 various constants 










si = , S2 = -1 , ri = 
ai = 1 , a2 = 1 , A = 


^4, r2-- 

1, P2 = 


= 1 
= 3 . 


Using 


(3.13) one 


obtains the following mass spectrum 










M2 = ^m(m + 3) + 0(m°) , 


m > 1 , 



(3.39) 



(3.40) 

where, it turns out, the constant ^ is still given by the corresponding expression in (3.19). In 
the limit when 61,62 ^ Po, the singularity structure of Eq. (3.10) changes [9]. In this limit 
however, the equation coincides with the equation for the 0"*"^ glueballs, and therefore the 
mass formula should be changed to 

M^ = — (m + l)(m + 3) + C>(m°) , m>l, (3.41) 

corresponding to the mass formula (3.19) of the O"*""*" glueballs shifted by one (since the lowest 
state should correspond to the zero mode of Eq. (3.10) and not to a glueball state [17]). This 
will be the formula used for comparison to the numerical results for a ^ mq. 
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3.2 Numerical evaluation of the mass spectra 

In the following, we present the results of the numerical evaluation of the mass spectra 
corresponding to the states described in Section 3.1. For every state, we will illustrate the 
dependence of the masses on the angular momentum parameter along a generic direction 
(chosen to be ai = 2a2), along the special direction ai = 02, and a table comparing the 
numerical and WKB results (and the lattice results for glueball states). 

3.2.1 Masses of the 0++ glueballs 

The equation for the 0"^"^ glueballs can also be written as 



du 



l.i^(u''-at){u'-at)-u'4)dj{u) 



- k'u-'fiu) = 0. (3.42) 



This equation is symmetric under the interchange of ai and 02, and reproduces Eq. (2.14) of 
Ref. [8] for 02 — > 0. This differential equation can be solved numerically using the shooting 
method as described in Ref. [4]. We require that the solution be normalizable (that is for 
u ^ 00 f{u) should vanish), and regular at the horizon uh- These conditions restrict the 
possible values of M^ to a discrete set, which can be identified with the glueball masses. The 
analysis of [8] demonstrated that the O"*"^ glueball masses are very stable against the variation 
of a single angular momentum parameter. The numerical solutions of Eq. (3.42) show that 
this statement remains valid for the whole range of angular parameters (ai, 02), except in the 
region ai = a2 ^ uq. This is consistent with the fact that in the WKB approximation (to 
order 1/m) the ratio of masses are independent of ai, 02 everywhere except at ai = 02 ^ "^^o, 
where the approximation breaks down. As mentioned before, ai = 02, mq = is the special 
region where the inner horizon coincides with the outer horizon. As a result, the factor 
multiplying the second derivative term in (3.42) has a double zero (instead of simple zero), 
and the behavior of the solutions is different. 

In Fig. 3.1 we show the behavior of the lowest eigenvalue of Eq. (3.42). The valley along 
Oi = 02 is related to the fact that the differential equation (and the physics of the model) 
is symmetric under the interchange ai ^^ a2. Note that the function is smooth except at 
the point Oi = 02 = 00 (or oi = 02, Uq = 0). In Fig. 3.2 we show the behavior of the ratio 
of the glueball masses along the direction ai = 2a2, which illustrates the fact that along a 
generic direction (by a generic direction we mean that it does not asymptote to ai = 02) the 
glueball mass ratios behave just like for the case with only one angular momentum, that is 
they change only slightly and take on their asymptotic value very quickly. The fact that the 
behavior of the glueball masses does not change can be seen comparing Fig. 3.2 to Fig. (2.1) 
of Ref. [8]. Table 3.1 contains the comparison of the lattice results, the numerical solutions 
and the WKB results along a generic direction (chosen to be ai = 202). 

In Fig. 3.3 we show a ratio of masses along the special direction ai = a2, where the inner 
and outer horizons come together as ai/uo -^ 00 (this is also the region where the WKB 

14 




Figure 3.1: The unnormalized values of the 0"'"+ glueball mass (the lowest eigenvalues of 
Eq. (3.42)) as a function of the two angular momenta. Note that this function is smooth 
everywhere except in the region ai = a2 —^ oo. 

approximation breaks down). In this region the mass ratios behave very differently than 
anywhere else and depart from the lattice results (for example, the weak-coupling lattice 
value for Mq++*/Mq++ ioi N = 3 is about 1.74, which is notably bigger than the numbers of 
Fig. 3.3 at large a). 



3.2.2 Masses of the 0"+ glueballs 

The differential equation for ^ glueballs can be written as 



1 



w 



d^[-{u^ - at){u' - at)d^xe,iu)] = -M^ 



M 



al){u' 



u 



{u* - a\){u^ - ai) - ulu^ 



XeM . (3.43) 



The corresponding mass spectrum can be obtained using a similar numerical method as 
for the 0"*"^ glueballs. The dependence of the lightest 0""*" glueball mass on the angular 
momentum along a generic direction (chosen again to be ai = 202) is given in Fig. 3.4. One 
can see that while the masses are fairly stable against variations of the angular momentum, 
just like in the case of 02 = discussed in [8], the actual values of the mass ratios compared 
to 0"'"+ increase by a sizeable (~ 25%) value. The change is in the right direction as suggested 
by recent improved lattice simulations [20] . The actual asymptotic value of the mass ratio 
^°~+ = 1.59 is the same as for the 02 = case everywhere except very close to the region 
ai = 02 ^ Mo- Table 3.2 contains the comparison of the lattice results to the supergravity 
results evaluated using the numerical and the WKB methods. 

Just as for the case of the O"*"^ glueballs, this ratio behaves very differently along the 



special ai = 02 direction, and departs significantly from the lattice result 
as it can be seen in Fig. 3.5. 



'o+^ 



I lattice 



1.46 



15 



1.594 



1.592 



1.5 



1.586 




10 



Figure 3.2: The ratio of the 0"*"^* mass to the 0^^ mass along a generic direction, chosen 
here to be ai = 2a2 = a. Note, that the change in the ratio is tiny, and the asymptotic value 
of the ratio is the same as in Ref. [8] in the case of ai -^ oo, 02 = 0. 




50 100 150 200 250 300 

Figure 3.3: The behavior of the ratio r of the mass of the excited 0^^* glueball mass to 
the 0"'""'' mass along the line ai =02- Note, that along this direction the solutions behave 
very differently than anywhere else in the parameter space and depart significantly from the 
lattice results. 
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state 


lattice 


numerical 


WKB 


0++ 


1.61 ±0.15 


1.61 (input) 


1.55 


0++* 


2.8 


2.57 


2.53 


0++** 


— 


3.49 


3.46 


Q++*** 


— 


4.40 


4.37 


Q++**** 


— 


5.30 


5.28 


Q++***** 


— 


6.20 


6.18 



Table 3.1: The masses of the first few O"*"^ glueballs in GeV. The first column gives the 
available lattice results [18, 19, 20], the second the asymptotic value of the supergravity 
calculation using the numerical method (the point is chosen to be ai = 2a2 = 20-Uo)5 while 
the third column the WKB result for the same supergravity approximation. 
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1.55 
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1.45 










1.4 


1 










1.35 


/ 












/ 2 


4 


6 
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10 


1.25 
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Figure 3.4: The ratio of the lowest 0^^ mass to the lowest 0^+ mass along a generic direction, 
chosen here to be ai = 2a2 = a. Note, that the ratio is very stable against the variations of 
the parameters. The actual change in the ratio is sizeable, and independent of the direction 
chosen in the (01,02) parameter space (except the line oi = 02) and agrees with the ratio 
found in Ref. [8] for the case of ai -^ 00, 02 = 0. As explained in the text, this figure is 
only reliable for the regions a <^ mq and a^ Uq which are shown by solid lines, while for the 
intermediate region denoted by a dashed line there are corrections due to the non-vanishing 
off-diagonal components of the metric. 
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1.25 

1.225 

1.2 

1.175 

1.15 
1.125 

1.075 



10 



12 



Figure 3.5: The behavior of the ratio r of the mass of the lowest 0^+ glueball mass to the 
lowest 0^+ mass along the line Oi = 02- As explained in the text, this figure is only reliable 
for the regions a <^ Uq and a ^ Uq which are shown by solid lines. 



state 


lattice 


Num.(ai,2 = 


= 0) 


Num.(ai = 2a2 = 


= 20) 


WKB 


0-+ 


2.59 ±0.13 


2.00 


2.57 


2.53 


0-+* 


3.64 ±0.18 


2.98 


3.49 


3.46 


0"+** 


— 


3.91 


4.40 


4.37 


g 1-*** 


— 


4.83 


5.30 


5.28 


g 1-**** 


— 


5.74 


6.20 


6.18 


g h***** 


— 


6.64 


7.10 


7.09 



Table 3.2: The masses of the first few O^'*" glueballs in GeV. Unlike 0^+ glueballs, the su- 
pergravity masses for these glueballs are sensitive to the values of 01,02. Two cases are 
displayed, illustrating the typical values that one gets for small and large Oi, 02 (the asymp- 
totic values for large 01,02 are the same for any generic direction, i.e. with oi 7^ 02). The 
lattice results are from [19, 20]. 



3.2.3 Masses of the KK modes of S^ 

In terms of u = ^/p, equation (3.29) for the first KK doublet (3.5) reads 



du 



u 



[{u'-at){u'-al)-u\l)dj{u) 



k 



0. 



where 



M-) = (1-4(1-4-^, 



w* 



H{u) = 4:ho{u) (1 



w* 



+ 



w 



:i- 



^2\2 



(3.44) 



(3.45) 
(3.46) 



The components of the second doublet (3.6) give the same equation with ai and 02 in- 
terchanged. Finally, the equation that determines the mass spectrum of the singlet (3.7) 
is 



M^a, 



u ^ 



u 



fin) 



u 



k'u' - 16m^)/(m) = 0. (3.47) 



This is symmetric under the interchange of Oi and 02. One can again numerically determine 
the solutions of these equations using the shooting method. In Figs. 3.6 and 3.7 we show 
the behavior of the 5*0(2) x 5*0(2) singlet mode first along a generic direction (which was 
again chosen to be Oi = 202), and then along the special direction ai = 02. One can see that 
this mode does not decouple on any region of the parameter space. Figs. 3.8 and 3.9 show 
the similar plots for the non-singlet KK modes (for (3.44)), which similarly do not decouple 
anywhere in the parameter space. Tables 3.3 and 3.4 show the comparison of the first few 
KK modes evaluated using the numerical and the WKB methods. 

3.2.4 Masses of the KK modes on the circle 

Next we consider the KK modes coming from the compact D-brane coordinate. These modes 
have the form (3.4), where (f){u) obeys the differential equation 



<9„ 



- ({u' - at)iu' - at) 
u ^ 



u 



du<Piu) 



[UjU 



e 



A^{u^ - u\ai + ai) - u^ + a^ai) 

(3.48) 



One can again numerically solve these equations. For a generic direction (chosen to be again 
ai = 2a2) we find that these modes decouple very quickly from the spectrum, just like in 
the case with one angular momentum parameter discussed in Ref. [8]. This is illustrated 
in Fig. 3.10. For the case of the special direction ai = a2, the numerical analysis of the 
decoupling is inconclusive. The masses of these KK modes grow much slower than for the 
generic case. At the point when our numerical solutions become unreliable, these modes 
are not decoupled yet, however one can not rule out the possibility that for a — » 00 they 
eventually do decouple (see Fig. 3.11). 
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Figure 3.6: The mass ratio r of the 5*0(2) x 5*0(2) singlet KK mode to the lowest 
glueball along the generic direction Oi = 2a2- 



1.4- 



1.3- 



1.2 




Figure 3.7: The mass ratio r of the 5*0(2) x 5*0(2) singlet KK mode to the lowest glueball 
mass O"*"^ along the special direction ai = 02. 
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state 


numerical 


WKB 


KK 


2.15 


2.00 


KK* 


3.23 


3.09 


KK** 


4.20 


4.09 


KK*** 


5.15 


5.05 


XX**** 


6.07 


5.99 


KK***** 


6.99 


6.91 



Table 3.3: The masses of the first few singlet KK modes in GeV. The first column gives 
the asymptotic value of the supergravity calculation using the numerical method (the point 
is chosen to be ai = 2a2 = 20mo); while the second column the WKB result for the same 
supergravity approximation. 
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Figure 3.8: The mass ratio r of the 5*0(2) x 5*0(2) doublet KK mode to the lowest glueball 
masses along the generic direction ai = 2a2- 
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Figure 3.9: The mass ratio r of the hghtest glueball to the 5*0(2) x 5*0(2) doublet KK modes 
along the special direction oi = 02- 



state 


numerical 


WKB 


KK 


2.84 


2.19 


KK* 


3.80 


3.34 


KK** 


4.73 


4.37 


KK*** 


5.54 


5.36 


XX**** 


6.57 


6.31 


KK***** 


7.47 


7.25 



Table 3.4: The masses of the first few non-singlet KK modes in GeV. The first column gives 
the asymptotic value of the supergravity calculation using the numerical method (the point 
is chosen to be ai = 2a2 = ^Ouq), while the second column the WKB result for the same 
supergravity approximation. 
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0.5 1 1.5 2 2.5 

Figure 3.10: The mass ratio r of the hghtest glueball to the KK mode on the compact D- 
brane coordinate 62 along the generic direction ai = 2a2 = a. Just as for the case with only 
one angular momentum, these states decouple very quickly from the spectrum. 




Figure 3.11: The mass ratio r of the lightest glueball to the KK mode on the compact 
D-brane coordinate 62 along the special direction ai = 02- 
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state 


numerical 


WKB 


KK 


11.27 


11.24 


KK* 


11.48 


11.45 


KK** 


11.76 


11.72 


KK*** 


12.09 


12.06 


j^j^**** 


12.48 


12.45 


KK***** 


12.92 


12.89 



Table 3.5: The masses of the first few KK modes on the circle in GeV. Along a generic 
direction these KK modes decouple, thus we have chosen an arbitrary point ai = 2a2 = 2.5mo 
for the comparison of the numerical and WKB results, which are given in the first and second 
column. 

4 Conclusions 

In this paper we have presented a two-parameter family of supergravity models of non- 
sup ersymmetric 3 + 1 dimensional SU{N) Yang-Mills theory, based on regular geometries 
with D4-brane charge. In these models, we have evaluated the mass spectra of the scalar 
glueballs and some of the related KK modes everywhere on the two dimensional parameter 
space using both numerical and analytic (WKB) methods. We find that the glueball mass 
ratios are very stable against the variation of the angular momentum parameters. The 
asymptotic values of these ratios for large angular momenta are in good agreement with 
the most recent lattice results everywhere in the parameter space except along a special 
line Oi = a2 3> Uq (which is exactly the region where the WKB approximation breaks 
down, and also the region where the inner and outer horizons approach each other and 
the supergravity approximation approaches a discontinuous limit). The KK modes on the 
compact D-brane coordinate decouple for large angular momenta everywhere (except perhaps 
along ai = 02, where our analysis of decoupling is inconclusive). The KK modes on the 
S'^ however do not decouple from the spectrum anywhere in the parameter space in the 
supergravity approximation used in this paper. 

The masses evaluated in this paper in the supergravity approximation can in principle get 
large corrections when extrapolating from the strong coupling (large A) regime to the weak- 
coupling regime of the Yang-Mills theory. If the spectrum of the corresponding string models 
at small A do indeed reproduce the Yang-Mills spectrum, a natural question to ask is why the 
glueball masses (or perhaps only the glueball mass ratios) would get small corrections, while 
the KK masses get large corrections. Since in the limit 01,02 ^ Uq the metric approaches 
the supersymmetric space (2.31), it is possible that a subset of the masses may be protected 
by supersymmetry (as explained in Sec. 2.3, the naive limit uq ^ does not lead to a 
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super symmetric theory, since it does not have fermions in the spectrum). A problem of 
interest is thus to investigate the supersymmetric model with uq = 0, and determine which 
scalars belong to short BPS multiplets, and which ones are in long multiplets. Since the 
masses of the scalars belonging to short multiplets should not be changed in the A <^ 1 
regime, this could explain why the ^ glueball masses are so close to the lattice values, 
and it may be used as a highly non-trivial quantitative test of the conjectured relation of 
supergravity to non-sup ersymmetric SU{N) Yang-Mills theory. 
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